Time scales in nuclear giant resonances 
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We propose a general approach to characterise fluctuations of measured cross sections of nuclear 
giant resonances. Simulated cross sections are obtained from a particular, yet representative self- 
energy which contains all information about fragmentations. Using a wavelet analysis, we demon- 
strate the extraction of time scales of cascading decays into configurations of different complexity of 
the resonance. We argue that the spreading widths of collective excitations in nuclei are determined 
by the number of fragmentations as seen in the power spectrum. An analytic treatment of the 
wavelet analysis using a Fourier expansion of the cross section confirms this principle. A simple rule 
for the relative life times of states associated with hierarchies of different complexity is given. 
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Nuclear Giant Resonances (GR) have been the subject 
of numerous investigations over several decades [l[ . Some 
of the basic features such as centroids and collectivity (in 
terms of the sum rules) are reasonably well understood 
within microscopic models @,||. However, the question 
how a collective mode like the GR disseminates its energy 
is one of the central issues in nuclear structure physics. 

According to accepted wisdom, GRs are essentially ex- 
cited by an external field being a one-body interaction. 
It is natural to describe these states as collective lp-lh 
states. Once excited, the GR disseminates its energy via 
direct particle emission and by coupling to more com- 
plicated configurations (2p-2h, 3p-3h, etc). The former 
mechanism gives rise to an escape width, while the latter 
yields spreading widths (r^). An understanding of life- 
time characteristics associated with the cascade of cou- 
plings and scales of fragmentations arising from this cou- 
pling (cf 0, remains a challenge. A recent high- 
resolution experiment of the Isoscalar Giant Quadrupole 
Resonance (QR) [1, H, Qi| provides new insight for this 
problem. 

It has been shown by Shevchenko et al. @ that the 
fine structure of the QR observed in (p,p') experiments is 
largely probe independent. Furthermore, a study of the 
fine structure using wavelet analysis [III E2, EH reveals 
energy scales [j| E3 in the widths of the fine structure 
displaying a seemingly schematic pattern, as can be seen 
in Fig. 2. This pattern varies with the structure of the 
nucleus being studied. While the physical meaning of 
the results of such an analysis is still being debated, we 
try here to offer a general explanation. However, we do 
not embark on a specific microscopic analysis, but rather 
make use of general and well-established techniques of 



many-body theory. Gross effects due to nuclear deforma- 
tion and coupling to the continuum Q are not discussed; 
we rather focus on the decay of the QR into configura- 
tions of various complexity. 

To proceed we use the Green's function approach. A 
central role is played by the self-energy whose finer struc- 
ture is imparted upon the Green's function via the solu- 
tion of Dyson's equation which reads [l4| 
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where we assume G°(lu) = S a .p/(uj — e) to be diagonal in 
the basis a, (3, . . . while the complicated pole structure of 
G(ui) is generated by that of the self-energy T, a ^(u). The 
pole structure of G carries over to the scattering matrix 
given by 

T a ,p(u) = £ Q ,,3M + Sa^'MG^a'MlV^M (2) 

from which a cross section ~ |T Q . Q (a;)| 2 is obtained. 

Within the excitation energy range of the QR the nu- 
cleus has a high density of complicated states of sev- 
eral tens of thousands per MeV and even more for heavy 
nuclei. These many states appear in the self-energy as 
poles in the complex energy plane close to the real axis. 
The small widths imply they are long-lived states and 
traditionally classed as compound states. The simpler 
intermediate structure of the excitation is expressed by 
the substantial fluctuations of the corresponding residues 
associated with the poles of the self-energy T,(ui) [lol ]. 
In other words, while the individual pole positions of 
S(w) are virtually unstructured (l6j . it is the variation of 
the corresponding residues that bears all the information 
about intermediate structure. Note that our approach 
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FIG. 1: Residues of the self-energy: for 6 (top) and 11 (bot- 
tom) intermediate states. The randomisation is clearly dis- 
cernible. 



differs from a traditional microscopic calculation in that 
we start from the outset from a random distribution of 
pole terms representing compound states. Traditional 
microscopic approaches cannot do justice to such struc- 
ture [ItJ and usually suffer from necessary truncation of 
configuration space and, associated with it, from possible 
inconsistencies and spurious states. 

We assume that the QR being a collective lp-lh state 
decays via a cascade progressing through (2p-2h)-, (3p- 
3h)-configurations and so forth to the eventual compound 
states. In turn, each of the intermediate states (including 
the initial QR) can either decay directly to the ground 
state or via some more complicated intermediate state. 
Below we will show that it is this mixture that is seen in 
the cross section and extracted by wavelet analysis, and 
it is the variety and cascading complexity of states that 
invokes the structure of the residues of the poles of the 
self-energy. Of importance to note is that the number 
of states available within the energy domain of the QR 
increases with its complexity: for example, six (2p-2h)- 
states, eleven (3p-3h)-states, down to several thousand 
compound states (the numbers six or eleven should be 
taken as examples without claim for quantitative cor- 
rectness). Moreover, the corresponding life times are ex- 
pected to increase in line with their increasing complex- 
ity, which is in accordance with their decreasing spread- 
ing widths (below we come back to this particular aspect 
of scaling) . 

As a typical case study we investigate here a wavelet 
analysis of a simulated cross section that results from a 
particular input for the self-energy. Since arbitrary units 
are used, we concentrate on the energy interval [0,1] and 
use for the pole position e = 0.5 — i0.5 of the single pole 
of Go |U). The number of compound states is assumed 
to be 300; this is of course much less than the experi- 
mental level density in the region of a QR for a medium 
or heavy nucleus, but it suffices for our demonstration. 
The real parts of the pole positions are assumed to be 
randomly distributed with a uniform distribution of the 
mean distance 1/300; the imaginary parts are randomly 
distributed in the interval [0.004,0.007]. 

For the sake of illustration we consider four sets of 



residues 

4 fi 2 

with an overall strength s = 10~ 5 . This order of mag- 
nitude is based on the mean value of the widths of the 
compound states being about 10~ 4 to 10 -5 times smaller 
than the T* With these residues the self-energy 

reads 

300 

EM = V (4) 

The poles at the complex positions ojk occur in the lower 
w-plane with to being the energy variable; the other sym- 
bols are explained in the text. If only i = 1 was to occur 
with /i = 6, a typical pattern of the residues = hi t k 
would be illustrated by the top of Fig.l; similarly for 
f-2 = 11 by the bottom. The inclusion of further terms 
would simply add additional peaks to the pattern. In the 
case presented below we have chosen fi = 11, fs = 17 and 
fi = 29 totalling to 6+11+17+29 additional peaks (not 
easily visualised, but beautifully discernible in the final 
analysis). We stress again that the four values fi were 
chosen for demonstration, more than four or other values 
can be used just as well. 

These arbitrary numbers used in the example chosen 
describe particular fragmentations of the QR into alto- 
gether 6,11,17 and 29 states of increasing complexity. 
The widths 7, giving rise to the Lorentzian shape of the 
residues are in reality determined by the product of the 
density of the compound states and the coupling of the 
i-th group to the compound states. The widths are the 
spreading widths of the respective states considered [l5| . 
As the complexity increases with label i we shall assume 
71 > 72 > 73 > 74- In the simulation we endow each 7$ 
with a random fluctuation with mean value 7i/4. As 
stated above we refrain from specifying a microscopic 
structure causing the residue pattern assumed for the 
self-energy; below it becomes clear that guidance comes 
from experiment. 




FIG. 2: Simulated cross section (left) and power spectrum 
(right). Units are arbitrary. The abscissa of the cross section 
is the unit energy interval, the energy values S on the abscissa 
of the power spectrum refer to the wavelet parameter using 
the same energy units. 

We also assume that each set fi uniformly distributed 
over the whole energy interval. This is similar in spirit 
to the assumption used in the local scaling dimension 



3 



0.03 
0.02 
O.Ol 




0.05 



O.l 



FIG. 3: Power spectrum for a particular asymmetric situation 
discussed in text. The dotted curve originates from a scan of 
the interval [0,0.5], the dashed curve from [0.5,1] and the solid 
curve from the total interval. Note that the most left peak is 
virtually absent in the dashed curve while fully present in the 
dotted curve. Units as in Fig. 2. 



approach [f|. The positions Pi in Eq.© are set to be 
~ 1/ fi which spreads the actual j ■ pi positions equidis- 
tantly over the whole interval with j running from 1 to 
however, we endow them with a small random fluctuation 
with mean value Pi/8. Note that the random fluctuation 
of widths and positions generate a mild degree of asym- 
metry in the energy interval [0,1], resulting in slightly 
different patterns in the intervals [0,0.5] and [0.5,1]. The 
near equality of the positions, that is - apart from slight 
random fluctuations - the regular pattern of the various 
fragments as illustrated in Fig.l, is basically dictated by 
experimental findings: if there is no near regular pattern 
there will be no discernible structure in the power spec- 
trum of the wavelet analysis. However, we shall return 
below to the case where such regular pattern may occur 
only in a smaller portion of the interval. 

The first obvious choice for the widths assumes simply 
7; = l/(2/j) yielding the simulated cross section shown 
in Fig. 2 (below a precise analytic expression confirming 
the l/(2/i)-law is given). A variation of such choice is 
rather significant, we shall return to this aspect in detail. 

The analysis using a Morlet-type mother wavelet (a 
contour plot of the wavelet analysis is illustrated in Fig. 5) 
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gives the power spectrum shown in Fig. 2; if not indicated 
otherwise we use the value k — 6 for the wave number of 
the mother wavelet. There is in fact a fc-dependence of 
the positions of the maxima of the power spectrum, it is 
given in analytic terms below. 

On the right part of Fig. 2 we clearly discern the four 
maxima that are produced by the four different values fi 
of the number of fragmentations. In fact, the fragmen- 
tation into fx = 6 produces (for k = 6) the maximum 
roughly at <5™ ax = 1/fx = 1/6; similarly, the other three 
maxima occur at 5f ax = 1//,:, i = 2,3,4. This is one of 
our major findings: 

the maxima of the power spectrum occur at 

sr x « fc/(27r) • 1 1 Si 

with I being the interval of the whole range of the QR 



considered and fi the number of fragmentations. The fac- 
tor k/(2it) originates from the analytic expression given 
in (7) below. 

The asymmetry found in some experimental data can 
obviously be accounted for by our analysis. We refer to 
cases where the analysis yields a pattern in the first half 
of the whole resonance being different from that in the 
second half, or in principle for any subdivision of the 
whole resonance. For illustration, we take fi = 14 while 
leaving all other parameters unchanged. In this way the 
total of 29 maxima of the residues r/ 4 are confined to 
only 14 within the left half of the interval. The effects 
are clearly seen in Fig. 3. Note that the positions of the 
maxima still remain unchanged. This type of asymmetry 
is clearly discernible in Fig. 9 of RefQjJ from the two- 
dimensional wavelet transform the wavelet power would 
give a similarly different pattern when taken at different 
portions of the whole interval. 

The folding (integration) of the cross section with the 
Morlet wavelet has to be done numerically. In order to 
obtain an analytic expression relating the number of frag- 
mentations fi to the positions of the maxima of the power 
spectrum, we consider an expansion of a cross section into 
a Fourier series 
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with the bulk term 
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(further terms with cos(mnuj / 1) are immaterial for the 
discussion). An intermediate structure manifests itself, 
if a few terms in ^ are appreciably stronger than the 
others. In Fig. 2 the terms with en ~ C22 ~ C34 ~ C58 are 
dominant; of course, terms for different m-values also oc- 
cur but are smaller by roughly an order of magnitude or 
more (here our analysis does not focus on m < 4: while 
giving larger contributions such values would correspond 
to 5 > 0.5 and represent gross and bulk structure). Per- 
forming analytically the wavelet-transform of each term 
in ((BJ (Mathematica gives a closed expression for the in- 
tegral from which the formula below can be extracted), 
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FIG. 4: Power spectrum: dependence of height at the maxi- 
mum on spreading width. The curve with the lower value of 
the most left maximum is identical to the one on the r.h.s. 
in Fig. 2., while the higher peak is due to a decrease of its 
spreading width or an increase of its life time. Units as in 
Fig.2. 
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one obtains an analytic evaluation of the positions and 
heights of the maxima of the power spectrum. For each 
sin(TO7ra;)-term the positions of the local maxima in the 
power spectrum turn out to be 



Max„ 



- + P 
2m7r 



I. 



(7) 



For k = 6 (and the unity interval I) this yields 0.16, 
0.088, 0.057 and 0.033 for m = 12, 22, 34 and 58, respec- 
tively as verified in Fig. 2. Note that a different choice of 
k moves the positions of the local maxima, yet the ~ 1 /m 
law prevails. The expression Q provides an obvious tool 
to be used to ascertain the number of fragmentations 
when the maxima are determined from an analysis of 
experimental data. Clearly, the number f m of fragmen- 
tations introduced above is related to the value m in © 
by m = 2f m . 

Furthermore, an increased value of k can resolve a peak 
in the power spectrum that is caused by two near values 
of /j. In fact, the distance between adjacent maxima (say 
rn = 17 and m — 18) roughly doubles when k is doubled. 

While - for fixed k - the 1/fi dependence of the max- 
ima of the power spectrum is an important finding, even 
more significant is the result that the values at the max- 
ima (the heights) also obey the same l//j-law if the cor- 
rcsponding Fourier coefficients are about equal. Indeed, 
a straight line can be drawn through the maxima in Fig. 2 
as the four values c m , m = 12, 22, 34, 58 are about equal. 
We recall that, for example, sin(127rx) generates //. = 6 
peaks of a width 7^ = 1/(2/^) in the energy (unit) in- 
terval for the cross section. This can be exploited in a 
realistic analysis: a deviation from this straight-line-rule 
signals effectively a deviation from the spreading width 
being assumed to be l/(2/j). This is illustrated in Fig. 4 
where the spreading width 1/(2/4) has been decreased 
to 1/(2.8/4). As a result, the value of the first peak be- 
comes enhanced. Since the spreading width is related to 
the life time of the states, we conclude: the life times are 
proportional to fi if the heights of the maxima lie on a 
straight line; an increased (decreased) height signals an 
even longer (shorter) life time. 

In this context we note that the number of peaks and 
troughs in Fig. 5 on the horizontal lines matches exactly 
the values of the /,: six on the top, further down eleven, 
then seventeen and twenty nine on the bottom. The 
actual values of these peaks and troughs determine the 
heights of the bumps in the power spectrum, that is the 
information about the life times of the respective frag- 
mented states. A similar wavelet transform obtained 
from experimental data is presented in Fig. 8 (and 9) in 
Ref|l0|; note that our schematic 'in vitro' illustration is 
of course much more symmetric. 

While in experiments the chaotic nature of the nucleus 
usually shows at higher excitation energies [l6j], the per- 
tinent structure revealed in the analysis may come as a 
surprise. We are of course familiar with order in the nu- 
clear many body system as shown in shell effects and 
simple collective states. The fragmentations of the QR 
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FIG. 5: Wavelet contour plot of cross section shown in Fig. 2. 
The symbols 5 and uio refer to the Morlet wavelet parame- 
ters used in ([5]). The (positive) maxima are in light shading 
and the (negative) minima in dark. For the top pattern the 
contours range from 0.4 to —0.4. 



may be due to a different quality: it could be a man- 
ifestation of self-organising structures [l8l . ITjjl [20l |. In- 
deed, the life-time of increasingly complex configurations 
of the QR is increasing toward the compound states and 
the ground state. There is no general accepted defini- 
tion of conditions under which the self-organising struc- 
tures are expected to arise. We may speculate that in 
the case considered here, once the nuclear QR state is 
created, it is driven to an unstable hierarchy of config- 
urations (metastable states) by quantum selection rules 
which connect these different complex configurations due 
to internal mixing. This problem needs of course a ded- 
icated study on its own and is beyond the scope of the 
present paper. 

We summarise the major points of our findings: (i) the 
position of the peaks in the power spectrum indicate the 
number of fragmentations of a particular intermediate 
state; the more complex states lie to the left of the simpler 
states (see Eq.(|7J)); (ii) the resolution of poorly resolved 
peaks can be improved by a higher value of fc; (iii) the 
values (heights) at the peaks are related to the spreading 
widths, implying knowledge about the life times: if they 
lie on a straight line, the life times are proportional to 
the number of fragmentations, if they lie above (below) 
the straight line the corresponding life times are longer 
(shorter). Finally, we mention that a pronounced gross 
structure of the experimental cross section as found in 
lighter nuclei, would have no effect upon our findings. In 
fact, such gross structure had to occur at the far right 
end (values of 5 appreciably larger than those used in 
the literature) of the power spectrum. 
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